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ABSTRACT 


In  this  paper  the  fracture  mechanics  of  orthotropic  materials  containing  collinear 
interface  cracks  is  considered.  The  primary  objective  is  to  study  the  influence  of  the 
thickness  and  the  structure  of  the  interfacial  regions  on  the  crack  driving  force.  The 
interfacial  region  is  assumed  to  be  a  relatively  thin  orthotropic  elastic  layer.  The  stress 
intensity  factors  or  the  strain  energy  release  rates  are  assumed  to  be  the  main  measure 
of  the  crack  driving  force.  .A  relatively  simple  and  efficient  method  is  presented  to  solve 
the  related  elasticity  problem.  The  results  are  obtained  for  a  wide  range  of  actual 
material  combinations.  In  order  to  study  the  influence  of  the  structure  of  the  interfacial 
zone,  the  problem  is  also  solved  for  isotropic  and  orthotropic  materials  bonded  through 
a  layer  with  hypothetically  selected  material  properties.  The  results  show  that  the  effect 
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of  the  thickness,  the  mechanical  properties  and  the  material  orientation  of  the 
interfacial  zone  on  the  strain  energy  release  rate  could  be  very  significant.  An 
interesting  and  a  rather  useful  result  obtained  from  the  collinear  crack  solutions  is  that 
the  strain  energy  release  rates  for  multiple  cracks  in  bonded  orthotropic  materials  with 
or  without  an  interfacial  layer  may  predicted  by  using  the  results  obtained  for  an 
isotropic  homogeneous  plane  provided  the  normalization  factors  are  selected  properlv. 

1.  Introduction 

In  recent  past  concerns  with  the  mechanical  failures  initiating  largely  at  the 
interlacial  regions  in  such  multiphase  materials  as  modern  composites,  thermal  barrier 
coatings  and  a  variety  of  other  bonded  materials  have  led  to  rather  extensive  studies  for 
the  purpose  of  understanding  the  interaction  between  the  flaws  that  may  exist  in  these 
regions  and  the  applied  loads  and  other  environmental  factors  [1],  [2]  and  [3],  From  the 
viewpoint  of  applications  generally  the  optimal  design  of  interfacial  regions  involves 
tradeoffs  between  strength  and  toughness.  For  example,  it  is  known  that  the  untreated 
carbon  has  very  poor  adhesion  to  epoxy  and  good  adhesion  can  be  accomplished  by 
treating  its  surface  and  by  using  proper  coupling  agents  (e.g..  variety  of  silanes)  [4]. 
However,  joints  of  such  high  interfacial  strength  have  usually  rather  poor  toughness. 
.\lso  ccuisiderations  regarding  load  transfer  and  fiber  integrity  require  a  relatively  weak 
or  compliant  interface.  Thus,  m  studying  the  failure  process  of  the  interfacial  regions 
the  mechanical  and  strength  parameters  are  expected  to  play  a  major  role.  In 
particular,  in  fracture  related  failures  the  crack  driving  force  would  heavily  depend  on 
the  constitutive  properties  and  the  thickness  of  the  interfacial  zone,  the  size,  location 
and  the  orientation  of  the  crack  and  on  the  properties  of  the  adherents  as  well  as  on  the 
iiiasnitude  and  the  nature  of  the  applied  loads  and  the  geometry  of  the  medium  [l]-[3]. 

Since  it  i‘;  not  likely  to  develop  a  process  that  would  produce  an  interfacial  zone 
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which  would  have  all  such  desirable  properties  as  high  strength,  high  toughness,  high 
fatigue  and  impact  resistance  and  high  resistance  to  moisture  penetration,  and  since  al’ 
these  properties  are  quantitatively  expressed  in  terms  of  continuum  concepts,  the 
fundamental  studies  relating  to  the  chemistry  of  adhesion  and  material  processing  must, 
therefore,  be  guided  by  the  corresponding  mechanics  considerations.  It  is  through 
studying  the  appropriate  mechanics  particularly  the  fracture  mechanics  of  the 
interfacial  zones  that  the  relative  importance  of  certain  material  properties  and 
dimensional  parameters  can  be  determined. 

In  most  studies  relating  to  the  fracture  mechanics  of  bonded  materials  it  is 
generally  assumed  that  the  composite  medium  is  piecewise  homogeneous.  On  the  other 
hand  a  close  examination  of  these  materials  seems  to  indicate  that  in  nearly  all  cases 
the  interfacial  region  has  a  structure  which  is  different  than  that  of  the  adherents.  In 
some  cases  this  region  may  simply  consist  of  the  reaction  zone  and  may  have  a 
thickness  of  only  a  few  lattice  parameters.  In  other  bonded  materials  the  interfacial 
regions  having  their  own  structure  and  thickness  may  develop  as  a  result  of  surface 
preparations  and  coupling  agents  used  in  processing.  For  example,  in  many  of  the 
j)olymer  matrix  composites  very  near  the  interface  crystallization  of  the  polymer  or 
cop(jlymerization  of  the  matrix  and  the  coupling  agent  seems  to  produce  a  highly 
oriented  region  with  a  columnar  or  lamellar  structure  which  may  generally  be  modeled 
as  a  thin  orthotropic  layer  (e.g..  PEEK-Carbon  <'omposites)  [4].  [9].  .A.mong  other 
examples  for  such  interfacial  regions  in  lioiuled  materials  one  may  mention  the  grain 
biuindaries  and.  at  a  larger  length  scale,  certain  geological  materials  such  as  shale  - 
sandstone  interfaces. 

The  main  objective  of  this  paper  is  to  study  the  influence  of  the  structure  and 
the  thickness  of  the  interfacial  regions  on  the  fracture  behavior  of  bonded  orthotropic 


materiab.  Specificallv.  the  effect  of  twM  tactors.  namely,  the  material  orthotropy  ami 
the  multijtle  sitt*  eraekius;  on  the  strain  eneray  release  rate  will  he  iin'estiaateil.  TliU'. 
the  tei  hni<[ue  iie\’eiope(l  is  suitahle  to  '-tmly  tht'  aeneral  prohlem  ot  i  ihline.ir  iutcrt.icial 
('rai'ks  in  honihal  urthotropii'  la\'ers.  In  this  papt'r.  howe'\'er.  it  will  he  assnmecl  th<it  the 
thickness  of  the  ailherents  in  comparison  with  the  crack  size>  ami  the  intertacial  zone 
thi<'kness  are  >ufficiently  larae'  so  that  the  perturhation  prolhem  considered  may  he 
approximateil  h\'  two  ilis>imilar  orthotropic  half  spaces  homh'd  thronsh  an  orthotrcjidc 
laver.  The  interface  crack  iirohlem  for  bonded  -emi-mfiiiite  orthotropic  plane-  was 
preidon-h'  cotisiilereil  in.  for  example.  :l()d,Tdh  .A.mona  ■'ome  of  the  siaiiificant  recent 
-tmlie-  (li'.dinu  v.’ith  the  interfact'  crai  k  prohh'ms  in  Itomleil  aiii-otropic  half  planes  one 
m<i\'  mi'iition  1-h-  lo  .  The  prohlem  ot  cracks  [lerpendicular  to  and  intersecting  the 
intt'i'face-  in  homled  orthotropic  lay«'rs  may  I)e  t<amd  in  ilC;. 


2.  ( )n  the  Formnlatitiii  ot'  th<'  Prohlem 

The  plane  -train  <  oUimvir  interface  <'rat  k  pridtlem  it)!'  two  orthotropic  half  planes 
1  and  3  homled  rhiouuh  an  orthotropic  iiiterfacial  layer  2  of  thickness  h  is  described  in 
Fiu.  1.  To  'implif\‘  the  formulation  of  the  problem  the  en^ineerins;  elastic  constants  E,j. 

and  G,^.  d.  i  =  1.  2)  for  an  orthotropic  plane  will  l>e  replaced  by  the  "averaged 
ela.-tic  parameter-  ilefined  by  17;.  'IS] 
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for  the  case  of  plane  strain.  Thus,  by  using  the  contraction 


yl=Xl/^^c,  y.^  =  x.^>lc. 


(3) 


and  defining 

’•n(yi-  y2)  =4- ^2)-  'I2(yi- >'2)  =  <^12(X1.  x^j).  -^(yi.  y.,)=C(7^^(Xi.  xj.  (4) 

'■i(yi-  y^)  =  u,(Xi.  xj.  V2(y,.  y.^)  =  u^lx,.  xj.  (5) 


the  stress  displacement  relations  may  be  expressed  as 

1  E. 
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dv, 


/  “ ■ 2  "'1 
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•2(^■o  +  :/o)  ^  dy-i  d}\’- 


(6) 


Here  x,.  .  and  u,.  (i,  j  =  1,  2),  are,  respectively  the  physical  coordinates,  stresses  and 

displacements.  By  substituting  from  (6)  into  the  equilibrium  conditions 


d'-ii  ,  d(Tj2  _  (, 
dX]  dXj 


d(T2i  &>T22  _  A 

dxi  ^x^  " 


(-) 


it  may  be  shown  that 


'dy.^Oy/  '^yi^y2 


0. 
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where 


2(«o  +  ‘"o)  j  _  2i^o(«o  +  '"o)  ,1 
,  2  •  '^  2  -  ;  2 


By  using  the  standard  Fourier  Transforms  the  solution  of  (8)  may  be  obtained  as 


follows: 


■dyi-yj)  =  ^  I  [^1^’^' ' 


+  C3e®'  ' ^  C^e"^^  I '  ''*]  da. 


'■2(yi-y2J  =  ^ 


-Si  I  a  y  >  .  .  -Sy  cv 


-^[  J3C,e-^i'“'>^  +  d,C,e 


^  .i3C3e®>  I ""  I  >■'!  -  I  “  ^  >"21  da. 


X  <  y,  <  cc  ,  0  <  >’2  <  h  ,  (10) 


where  the  functions  C,(a).  (i=1....4)  are  unknown. 


and  S]  and  S2  are  the  roots  of  the  characteristic  equation 


s'*  -  2«(jS^  +  1=0. 


Equation  (12)  has  four  roots  s,....  s^  given  by 


s,  =:  -s3  =  ^«:„  +  ^j«:,/-l  .  Re(s,)>0. 
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s  ,  =  -  S4  =  \  «o  “  \  ''o*  ~  1  )  >  0- 


(131 


Note  that  the  roots  s,...,  S4  depend  on  n-q  only  and  they  are  all  real  for  >  1  and 
complex  for  -1<Kq<  1.  For  «o  <  -  1  the  roots  are  pure  imaginary  for  which  equations 
(10)  imply  that  the  boundary  value  problem  has  no  feasible  solution.  This  restriction  on 
Kq  may  also  be  arrived  at  by  requiring  that  the  elasticity  matrix  C,j  (in  <t,  =  C.^tq)  be 
positive  definite  [19],  From  (1)  and  (2)  it  may  easily  be  shown  that  for  isotropic 
materials  s:,j  =  1  and  equations  (S)  and  (12)  reduce  to  the  known  results.  Formally. 
e(piations  (10)  give  the  solution  for  an  orthotn^pic  layer  where  the  constants  CJ....C4  are 
•  lefermined  from  the  stress  and/i^r  displacement  boundary  conditions  prescribed  at 
Xj  =  constant  planes. 

3.  Integral  Equations  for  Collinear  Cracks 

Consider  now  the  plane  strain  or  plane  stress  collinear  crack  problem  for  two 
orthotropic  half  spaces  bonded  through  an  orthotropic  layer  described  in  Fig.  1.  It  is 
assumed  that  the  elasticity  problem  for  the  composite  medium  has  been  solved  under 
given  applied  loads  by  ignoring  the  cracks  and  the  stresses  (7,2'’(x,.  0)  and  cr,2'^(Xi,  0) 
along  the  interface  X2  =  0  have  oeen  calculated.  Thus,  by  using  a  superposition  the 
i)riginal  crack  problem  can  be  reduced  to  a  perturbation  problem  in  which  the  crack 
surface  tractions  -it22°(Xi,  0)  and  -a-,2‘^(x,,  0)  are  the  only  external  loads.  Note  that 
this  process  is  rather  straightforward  if  the  stress  state  at  infinity  is  zero  and  the  applied 
loads  consist  of  body  forces  and  moments  having  finite  resultants  and  acting  within  a 
bounded  region  of  the  medium.  Otherwise,  the  stresses  acting  on  the  medium  at  infinity 
must  be  such  that  the  strains  along  the  interfaces  X2  =  0  and  Xj  =  -h  in  the  three 


materials  remain  <’umparible.  Fur  f’xample.  if  the  applie>l  Ktads  at  infinity  are  uiiifurm 


tension  =  (t,|  and  uniform  shear  to  maintain  |■umpati^iIiry  in  the  'tr.iiii' 

•  ,i  or  prevent  any  mismatidi  in  the  rh'^phieement  components  u,  on  x_,  =  (J  and  x,  =  -  h. 
it  Would  he  nei'essary  to  apply  to  material"  1  and  2  the  followimr  ad<liti(uial  external 


loaUs! 


I  ti 


-  V.  x,i  =  E,, 


,  F  '  ■ ' 

'  II  *-!  ! 


I  I  1 


'11 


E.,'*' 


—  I.  I)  X.  ■: 


1 J 


Eld-" 


—  I.  -h".  X, 


I  Ui 


when-  the  sujier'irijit'  i  1  i.  i2'  and  'oi  refer  to  mattu'ials  1.  2  and  3.  resperti\-(dy  iFin. 
1  I, 

The  formulation  oiven  in  the  previous  st'crioii  is  valid  for  each  of  the  three 
materials  shown  in  Fiu.  1.  An  indi'x  k  =  1.  2.  3  will  he  added  to  each  quantity  to 
desionate  rlu'  materials  1.  2  <ind  3.  re'-pectively.  That  is.  the  cpiantities  which  appear  in 
Section  2  will  he  riqdaced  hv  x^,,.  E^,,.  iq,„.  .q,.  iq.,.  cq,.  s^.,.  C^.j.  (h  =  1. 

2.  3:  I  -  1.  2:  /  =  1.  2;  j  =  1.  2.  3.  4i.  Referrinii  now  lo  i  10).  since  the  displacements  in 
the  medium  I'ani'h  at  x ,  -  ~  tjr  v ,  =  ~  we  must  have 


C,;  =  1).  C,,  =  ().  C,,  =0.  C,^  =  o. 


1 15 1 


The  remaiiiiini  einht  functions  .irv  then  derenmned  hv  usin«;  the  continuity  and 
houndaiu'  conditions  at  Xj  =  0  and  x,  =  -  h.  The^e  I'onilitions  may  lie  expressed  as  (Fi«;. 
ll 

-  I  .  ,i  X|  .  1 1  1  =  .T  ,  ,  ,1  X| .  I)  !.  1^1  I  ,1  X| .  0  i  =  iT  X| .  0  ).  -  V  <  X|  <  S".  !  16  I 
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~  li)  =  ^322(^1-  -h).  <^212(^1-  ~  h)  =  <T3,2(x,.  -h).  ->c<Xi<x;.  (17) 

U2i(x,.  -h)=U3,(x,.  -h).  U.,2(-X,,  -h)=U32(x,.  -h).  -X<Xi<>:.  (IS) 

<T,22(Xi-  0)  =  -0-22°(Xi,  0)  =  P,(Xi),  ^^,I2(Xl•  0)  =  -'^r2^(-'^l-  *^)  =  P2(^l  '■  -'<1  S  L.  (  19l 

Uii(Xi.  0)-n2i(x,.  0)  =  0.  Ui2(xi.  0)-u22(x,.  0)  =  0.  XiSL'.  (20) 

It 

L=Z  Lj-  =  (a^  <  Xi  <  b/).  j  =  1....  n.  (21) 

j  =  1 

where  pj  and  p,  are  known  fnrctions.  represents  the  jth  crack  (Fig.  1).  and  L'  is  the 
complement  of  L  in  (  -  x.  x). 

By  using  (lo)  and  substituting  from  (10)  into  the  homogeneous  equations 
(16) -(IS),  six  of  the  unknown  functions  may  be  eliminated,  the  two  remaining 
unknowns  are  then  determined  by  using  the  mixed  boundary  conditions  (19)  and  (20). 
To  reduce  ( 19)  and  (20)  to  a  system  of  integral  equations  we  define 

[  Uii(Xi.  0)  -U2i(x,,  0)]  =fi(X,).  -  X  <  X,  <  X. 

-^  [  Ui2(Xi.  0)-U22(Xi,  0)]  =f2(x,).  -X<x,  <X.  (22) 

From  (20)  and  (22)  it  follows  that 

fi.(x,)  =  0.  (k  =  1.2).  X,  eL'.  (23) 


9 


4(x,)  dxi  =0.  (k  =  L  2:  i  =  1,...  n). 


It  is  seen  that  by  substituting  from  (10)  into  (16).  (17).  (IS)  and  (22)  all  eight  unkinuvu 
functions  may  be  expressed  in  terms  of  f,  and  f_>.  Ecjuations  (19)  would  then  give  a 
pair  of  integral  equations  to  determine  f,  and  fj.  This  process  is  somewhat  lengthy  but 
straightforward.  Thus,  referring  to  [20]  for  analytical  details  we  obtain 


-,ifi(xi) +  *,,2  i  f  I  SI  K,  (x,.t)f^(t)dt  =  p,(xi).  X,  gL. 

*  j  _  I 

L  LJ 


■nh  +i  K2^(xi.r)f^(t)dt  =  p.2(Xi).  x,  G  L. 


L  J  = 


where  K,^.  (i.  j  =  1.  2)  are  known  bounded  functions  and  the  bimaterial  constants  7,^.  (i. 
j  =  1.  2)  are  given  by  (see  [20]) 


-  -11 
'ij  -  n,j 


(i.j  =  l.‘2). 


li  1  1  —  —  11  )  1  — 


i  "  ^10  ^  ‘ 


^2(l-^^-,„)  ^2(l+'':■2o) 


C2E.20 


C(\2(l  +  «'jo)  2(  1  + 


.  iIq  —  n,  ,1121  -  n,i 


The  elastic  constants  which  appear  in  (27)  are  defined  by  (1)  or  (2).  Note  that  712  and 
*•  ,,  are  alwa\‘,s  positive. 

The  integral  equations  {‘25)  may  further  be  simplified  by  defining 
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■■>v  1 


flX|l  =  ;;  't’jlx, 


pi  X,  i  =  l)jl  X,  )  -  (  j;  pji  X,  I.  X|  L. 
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TliU"  liy  u>iu2  I  2^  I  and  i  29  i  and  roiul'iniini  i  2oi  ir  may  easily  lie  >liu\vn  riiar 


1  1  ^  '  ’’  ' 

^  I  r  -  X, 


dr  -  fi  x. 


,K,I  X,.  r  I  t'l  r  I  -  K_.i  x,.  r  i  fir  i  dr 


=  -  pi  X,  I. 


x,eL. 


31)1 


wlieie 


K,ix,.  ri  dd,, -K,,  -Iii  K.,  -I/Ki.id 

I-'  ji 


K  a  X I .  r  I  =  -  ^  K , ,  ~  K , ,  -r  /I  ^  K , ,  -  ;/K , ,  id 

-\'  ij‘  n  ■  ■■ 
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Fur  rile  ease  uf  n  eullinear  eraeks  L,....  L„.  rheuretieally  rhe  index  of  the  singular 
inreiir.d  equarimi  i3()i  is  n  and.  eonsecjnenrly.  ir■^  'olnrion  contains  n  arbitrary  (complex) 
cun.'re.nr'.  Thi''e  ci,u>ranr>  may  be  olnained  from  the  ■'ineie-valuedness  conditions  124) 
v.dneii  i!ia\'  now  be  expressed  as 

/, 

/ 

fi  r  I  dr  =  I)  .  _/  =  1....  n.  (32) 

If 

/ 


4.  The  Merhi  id  of  Sulurjull 

Tirnii'i  ru  -olve  rhe  inreoial  eipniriuii  l3lli  which  is  defined  on  rhe  union  of  i'  arcs 
L|....  L„  direi  rl',-  i:  hiiihlv  impracrical.  It  is.  however,  rather  easy  to  reduct-  ( 30 1  to  a 

1 1 


system  of  n  simpler  integral  equations  each  of  which  is  defined  on  a  single  arc  Lj.  Thus, 
by  defining 

b ,  -  a ,  b ,  +  a , 

t  =  — -  s  +  i  .  aj  <  t  <  bj.  -  1  <  s  <  1  . 


bi.  —  at.  bt  +  ai.  ,  ,  . 

Xi  =  ^  2  ~  ^  ^  2  "  •  ^  <  bfc-  -l<r<l. 


(33) 


t(t)  =  a^  <  t  <  bj;  p(Xi)=qt(r).  a^  <  Xj  <  b^. 


(34) 


KdXi-  t)  =  R*,j(r.  s).  K2(xi.  t )  =  Sfcj(r.  s).  a*.  <  x,  <  b^..  aj<t<b^.  (35) 


l,  (  r  s  i  =  _ _  =  1 

■  (bj -aj)  s-lb^-a*)  r  +  b^+a^-bi. -a^  ^-^kj 


l36) 


from  130)  it  ma\'  easil}’  be  seen  that 


1 


1  .  n  ' 


^54^  ds  -  -  c^.ir)  ^  '  lH^(r.  s)e^(s)ds 


-  1 


J  =  1  1 


^  [Rfcj(r.  s)fj(s)  +  Sfc^(r.  s)i  j(s)]ds  =  iqfc(s). 

;  =  1  "  _  1 


k  =  1,..,  n.  -  1  <  r  <  1. 


(31 


where  the  prime  in  the  third  term  indicates  that  the  summation  excludes  the  term 
rorresponding  to  j  =  k.  Xot^  that  the  kernels  h<.j(r.  s)  as  well  as  R*.j  and  Si^j  are 
l)ouiided  ill  the  dosed  interval  -l<(r.  s)  <  1.  The  unknown  functions  cq.  and  known 
functions  ([f.  respectively  represent  the  crack  opening  displacements  and  the  crack 


surface  tractions  for  the  kth  crack  L*.. 


We  now  observe  that  the  weight  function  of  i  j.(s)  in  (37)  is  [20j,  [21 1 
W(r)  (1 -r)'^(l +  r)'^  .  -  1  <  r  <  1. 

where 


a  = 


<3Si 


( 39 ) 


The  orthogonal  polynomials  associated  %vith  w(r)  are  the  Jacobi  polynomials 
and  therefore  the  solution  of  (37)  may  be  expressed  as 


'■  i,(r)  =  gi.lr)w(r).  g*,.(r)  =  ^  I r i.  -l<r<l.  k  =  l....n.  (40) 

j  =  (i 

where  are  unknown  complex  constants.  From  (34)  it  may  be  seen  that  the 

conditions  (32)  become 


.1 

I.  r  )<.lr  =  0  .  t:  =  1....  n.  (41) 

_  I 

Thus,  by  usine:  the  orthogonality  conditions  for  the  Jacobi  polynomials,  from  (40)  and 
I  41 )  we  fill'] 

Afco  =  0.  k  =  l....  n.  (42) 

The  singular  integral  equations  (37)  may  be  regularized  by  using  the  following 
property  of  the  Jacobi  polynomials; 


1 


W(S)P,/^*"^'(S) 


ds 

s  —  r 


--  w(r)P„"'  -^’(r)  = 


-  1  <  r  <  1  . 


(431 
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Thus,  by  substituting  from  (43)  into  (37)  the  singular  terms  in  the  integral  equations 
may  be  removed  and  (37)  may  be  reduced  to  a  system  of  functional  equations  with 
(k  =  1....  n:  j  =  1,  2.  3...)  as  the  unknown  coefficients.  These  equations  may  be  solved  bv 
truncating  the  series  defined  in  (40)  at  j  =  N'n.  and  by  using  a  method  of  collcjcation.  Fur 
improved  accuracy  in  each  equation  the  collocation  points  r,„  are  selected  as 


T.vjrm)  =  0-  =  tf„  =  ^(2m-l).  m  =  ^  =  1....  n.  (441 

The  integral  equations  '37)  may  then  be  reduced  to  a  system  N'l + N'„  linear 
algebraic  equations  in  .4.^.^.  (k  =  1....  n.  j  =  L...  N\.).  It  should  be  noted  that  in  addition 
to  the  first  two  terms,  the  third  term  involving  the  kernels  h^.^  in  (37)  can  also  be 
evaluated  in  closed  form  [20].  Furthermore,  the  integrals  in  the  fourth  term  are  all 
Gaussian  type  and  may  be  calculated  rather  accurately  without  any  difficulty. 


5.  Stress  Intensitv  Factors  and  the  Strain  Enersv  Release  Rate 

It  may  be  seen  that  once  the  coefficients  A^.j  are  determined.  (40)  essentially 
provides  a  closed  form  solution.  From  a  view  point  of  fracture  mechanics  generally  the 
quantities  of  interest  are  the  crack  opening  displacement,  the  stress  intensity  factors  and 
the  strain  energy  release  rate.  From  (22).  (29).  (33)  and  (34)  the  crack  opening 


ilisplacement  may  be  expressed  as 


)  —  7j  (’-^12  ~  ’^22  )  +  '(^11 


~  l 

-U2,")=  f(t)dt 


U-  -  ar- 
o 


i,q.(s)ds.  k  =  1....  n. 


Substituting  now  from  (40)  and  evaluating  the  integral  [20]  we  find 
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I 


(rt  +  +  1 ) 


(r)  .  -  1  <  r  <  1, 


j  =  1 


r  =  ^-^ - k=l....n. 


2 


(46) 


In  the  mixed  boundary  value  problems  such  as  that  under  consideration  it  may 
be  shown  that  the  fundamental  solution  or  the  weight  function  w(r)  defined  by  (3S) 
represents  the  asymptotic  behavior  of  stresses  and  displacements  near  the  crack  tips.  In 
this  problem  w(r)  has  complex  singularities.  Thus,  one  may  define  a  complex  stress 
intensity  factor  as  follows: 


k(b,)  =  r, 

iki  +  ik2  =  /im  Nf2(x,-b*.) 

0)  +  ><^i 

.2(Xi.  0)]. 

k(afc)  =  >1 

ik,  +  ik2  =  ^  (afc  -  X, )  " 

0)  +  icr, 

2(X,.  0)1, 

k  =  l..., 

Observing  that  left  hand  side  of  (30)  represents  the  stresses  on  L'  as  well  as  on  L. 
referring  to  [20]  for  details  the  stress  intensity  factors  defined  by  (47)  may  expressed  as 
follows: 

=  -  72n  1  ~  ^  (b*.  -  a^)’“'gfc(l)^  (b^.  -  a*.)/2  , 


k(aO  =  721^1  -  7^  (bfc-aj  -  1)^(1)^  -  aO/2  . 


(48) 


where  gjt(r)  is  given  by  (40).  Similarly  in  [20]  it  was  shown  that  from  the  ‘‘crack  closure 
energy"  the  strain  energy  release  rate  may  l)e  evaluated  as 


g(b,)  = 


rk(br)k  (b^ 


•^1 


5(ajt)  = 


Tk(a^)k(ajt) 


'21 


15 


(49) 


where  the  elastic  constant  721  is  given  by  (26)  and  (27)  and  the  bars  refer  to  complex 
conjugates. 


6.  The  Special  Cases 

In  this  section  vve  give  the  results  for  some  simple  special  rases  for  which  the 
solution  may  be  obtained  in  closed  form. 


(a)  A  single  interface  crack  (  -  a.  a)  in  two  bonded  dissimilar  half  planes  (n  =  1. 
h  =  >c,  Fig.l) 

In  this  case  the  Fredholm  kernels  Kj  and  K2  in  (31)  are  zero  and  the  solution 
may  be  expressed  in  closed  form  [20].  In  particular  for  a  uniform  stress  state  at 

infinity,  or  for 


(722(Xi.  0)  =  -(Tg.  <T,2(X,.  0)  =  -Tg. 

vve  obtain 


(50) 


Ai  =  -7=J==(  7A==  +  i^).  Afc  =  0.  k  =  2.  3 . 


giving 


(51) 


(52) 


k(a)  =  (2a)''‘'(l  -  2iu.){r]crQ  + 


k(  -  a)  =  (2a)  “*'(1  +  2iw)(r;(7Q  +  tT-o)Aia. 


(53) 


nW  + 


(54) 


where  the  elastic  constants  are  defined  by  (26)-(28)  and  (39).  and  2a  is  the  crack  length. 
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If  the  medium  is  homogeneous  and  isotropic  with  the  elastic  constants  /i  and  k 
we  have 

-.=0.  -=0.  =  =  (j  =  1.2).  (59) 

where  is  the  shear  modulus  and  ft:  =  3-4e  for  plane  strain  and  «:  =  (3  - 1/)/(  1  + 1/)  for 
plane  stress  conditions  and  for  uniform  loading  kj  and  k^  are  given  by  (57). 

(c)  Bonded  isotropic  half  planes  with  a  single  crack  (  -  a.  a) 

If  the  two  half  planes  are  isotropic  having  the  elastic  constants  /i,.  >Cj  and  ^2. 
we  have 
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0  = 

4  4^.>  4/i, 


III  the  rase  of  uniform  applied  loads  (50)  the  stress  intensity  factors  are  again  given  by 
( 53 )  ivith  //  =  1. 

(d)  Homogeneous  isotropic  plane  with  two  collinear  cracks  of  equal  length 

To  provide  a  benchmark  for  the  collinear  interface  crack  results  we  give  the 
solution  for  an  isotropic  plane  having  two  collinear  cracks  of  equal  length.  Referring  to 
Fig.  1.  let  a,  =  -B.  bj  =  -  A.  a.2  =  A  and  b^  =  B.  For  uniform  tension  =  ctq  at 
x_,  =  X  the  mode  I  stress  intensity  factors  may  be  expressed  as  [22] 


,  [(BV.V)  E(k)/Kikl  -  1|  . 

ki(B)=  [1-E(k)/K(k)i  .  k*  =  l-(AVB')  .  (62) 

\|1  -  (A  /B  ) 


where  Kikl  and  E(k)  are  complete  elliptic  integrals  of  the  first  and  the  second  kind, 
respectively.  From  (62)  it  may  be  shown  that  the  limiting  values  of  k,(A)  and  k,(B)  are 


for  I  B  -  A  I  =  11  iii't .inr  B  -  A  '  —  v .  <tn<l 

ki(  B  )— .  k,iAi  —X  .  I  C4  • 

tor  B  =  I  > 'ii'r.iiir .  A— H.  In  rhr  l.irtcr  ut>rinu  rlmr  B  =  i  i 'ii>raut  and  koBi  i'  a 

tniic  rum  ct  A.  ir  m.i’''  aUn  ii>-  -ni.wn  fh.ir 

liiii  k,  B  =  -  X  , 

A-O 

Miaf  lx.  {,  iv  'luall  ','al’ii''  '4  A.  kj'B'  1'  luit  wt'll  •k'fiut'il.  Tin’  ia>irt'spi)Uilnis  'Train  onornv 
Ifli'a'f  r<i''''  aT  fill'  I  lai  k  !lli'  A  <inil  B  all'  UlVrn  liV 


Q  1  = 


k,dB 
1 1  - 


I 


\ 


Similarly  fur  i  B  -  A  x  we  h.ivu  rho  rofori’in’t'  value 


iG6) 


Si  =  g^  =  g* 


t^„-'iB  -  A  1/2 

S/I  M  1  -  e  I 


(67) 


, .  Xnmerii  ai  Rr'iilr-  and  Dixi  n>''ion 

Fur  rliu  liuiidi’d  dissimilar  marerials  with  inrvrface  rracks  from  (37)-(39)  it  is  seen 
rliat  rhe  inTcural  eiiuatioiis  are  uf  rlie  second  kind  ithat  is  -  91^01.  and.  consequently,  the 
funikiinental  sulutiim  tv  wliich  controls  asymptotic  nature  of  stresses  and  displacements 
near  the  crai  k  tips  lias  complex  siuoularities  resultins;  in  the  well-known  oscillatory 
-rri  S'  ,ind  di'placement  Wehavior.  For  example,  for  a  single  interface  crack  in  two 
liuiided  urr liutiupic  half  plane'  under  uniform  stre'xes  <it  infinity,  from  (52)  the  crack 
upeiiimi  di'[)lai  eiiieiiT'  may  he  uhraiued  as 


lit 


U,_;(X,.  +0l  -  Uj_,(Xi.  -Oi  = 


\ja--x,-  .  a, 3 


f—  I  <xjs(  _  log;  ( 


a  -  X, 


^0  '  1  H  —  Xi 

- MIIU  log  I  a^TYT 


\i  1-’  -ji 


II I  -  Uj, i  X,.  -  0 1  = 


\a- 


^1  '  'o  1  a  Xj 

==T-  I  X—  f'<)S(  *  log  I  ^  I  ) 


■  ji 


\  u'  -ji 


<iiil  -  log 


a  -  X, 
a  -  X 


l6S 


E<[uat;()ii  (6Si  ^how^  rhaf  as  Xj— Ta  the  ilispla«"eiiie:irs  wouM  oscillate  leailiiig  to  an 
.ippaieiit  crack  >urtai  e  iiitertereiice.  The  hiinaterial  coiisranr  _  which  controls  these 
ijM'illatiuiis  is  ilefiiieil  hy  i  39  I  in  t»'rius  <>t  -  .  which  in  turn  is  gieen  h}'  l26)-(2S)  or  (60l. 
It  i,'  perhaps  wiurhwhile  to  einplia>ize  that  in  crack  proldeins  there  are  really  two  kiinls 
ot  crack  surtace  intt'rtereucc.  One  is  the  mathematical  coiiseipience  of  the  assumed  jump 
discontinuities  in  material  constants  across  the  interface,  giving  oscillations  even  under 
pure  mode  I  loading  conditions.  However,  for  practical  material  pairs  the  size  of  this 
interference  region  having  the  forms  Au  ~  ^a  -  x,  •  sin(  -iog( a  -  x, ) )  or  Au  ~  .,ja  -  Xj 
■  cosi  _log|  a  -  X|  I  t  is  usually  of  the  order  »jr  smaller  than  the  microstructural  length 
parameters  and  its  effect  on  the  fracture  analysis  is  rather  negligible.  Furthermore,  this 
ItehavicU'  may  be  removeil  by  assuming  a  more  realistic  material  model  in  which  the 
jump  <lis<ontinuiries  are  replaced  by  "thin  interfacial  regions  with  steeply  varying 
materi.d  piajpertie-  j23j. 

The  second  kind  of  crack  surface  interference  is  macroscopic  and  geometric,  and 
is  invariably  the  consecpience  of  mode  II  loading  in  crack  problems  that  lack  symmetry 
with  respt'ct  to  tin'  plane  of  the  crack  (e.g..  an  embed<led  crack  parallel  to  the  interface 
or  the  free  -urfacei.  In  the  case  of  au  interface  cra<‘k  in  bonded  dissimilar  materials. 
'iiK  e  the  plane  of  the  crack  is  not  a  plane  of  symmetry,  under  pure  mode  II  loading  the 
normal  <  i»mi)oneut  of  the  crack  surfa<  e  opening  wouhl  be  an  o<ld  functii)n  resulting  in 

•JO 


surface  interference.  As  an  example.  Figure  2  shows  the  normal  component 
Av  =  Ui2  -  U22  ~  of  tlio  crack  opening  displacement  given  (68a)  for  two  bonded 
dissimilar  orthotropic  half  planes  for  various  relative  values  of  tensile  and  shear  loadings 
(Tq  and  The  calculations  neglect  the  interference.  Hence  along  a  substantial  portion  of 
the  crack  Av  seems  to  be  negative.  For  various  <Tg/r„  ratios  the  values  Xq  for  which 
Av(Xq)  =0  are  i^iven  as  follows: 

0  0.1  0.2  0.5  1.0 

.Xo/a:  0  -0.60536  -0.88271  -0.99693  -0.99996 

These  results  indicate  that  for  small  values  of  (Xy/r^  the  interference  is  macroscopic,  the 
solution  given  in  this  paper  is  not  valid  and  to  obtain  the  correct  solution  the  problem 
must  be  formulated  as  a  crack/contact  problem  (e.  g.,  [24]).  However,  these  results  edso 
show  that  if  the  mode  I  component  of  the  external  load  is  sufficiently  high,  the  possible 
crack  surface  interference  may  safely  be  neglected. 

From  (26)-(28)  and  (39)  it  may  be  seen  that  the  crack  tip  stress  and 
displacement  oscillations  disappear  if  7  =  0  or  if 

.  (69) 

b-jo  Fjo 

where  the  “averaged"  material  constants  E,q  and  e,o.  (i  =  1-  2)  are  defined  by  (1)  for 
plane  stress  and  (2)  for  plane  strain  conditions.  From  (69)  it  also  follows  that  in 
isotropic  incompressible  materials  Uq  =  1/2  for  plane  stress  and  eg  =  1  for  plane  strain 
and.  therefore,  oscillations  disappear  only  under  plane  strain  conditions. 


In  this  study  two  types  of  examples  are  given.  In  the  first  some  specific  materials 


are  considered  in  various  combinations.  The  elastic  properties  of  these  materials  are 
listed  in  Table  1.  Materials  a-c  are  fiber  reinforced  structural  composites  (graphite- 
epoxy).  Except  for  a  90  degree  rotation  about  the  z  axis,  materials  a  and  b  are  identical. 
Material  k  is  a  crystalline  solid  (Topaz)  and  again  material  /  is  the  same  as  k  except  for 
a  90  degree  rotation.  Materials  d-j  are  isotropic.  The  material  combinations  used  in  the 
examples  are  shown  in  Table  2.  In  this  table  materials  nii,..,  m4  are  assumed  to 
represent  the  grain  boundary  region  in  bonded  Topaz  crystals.  They  are  assumed  to  be 
isotropic  having  hypothetically  selected  elastic  constants. 

Table  .3  shows  the  material  combinations  used  in  the  second  group  of  examples  in 
which  two  dissimilar  isotropic  materials,  carbon  and  epoxy,  are  assun  ^’d  to  be  bonded 
through  an  orthotropic  interfacial  layer  of  hypothetically  selected  properties.  Depending 
on  the  processing  and  coupling  agents  used,  the  interfacial  zone  may  exhibit  either  a 
columnar  (E.222  >E2ii)  ^  lamellar  (E222  <E2u)  structure  (E222  ^211  being  the 

Young's  moduli  in  X2  and  x,  directions,  respectively).  In  each  example  only  one 
material  parameter  (E2n-  E222  or  ^212)  is  varied.  Since  three  of  the  Poisson’s  ratios  are 
fixed,  the  remaining  three  must  be  varied  along  with  the  Young's  moduli  to  maintain 
the  symmetry  of  the  elasticity  matrix.  The  table  also  shows  the  oscillation  parameter  -j 
defined  by  (39)  and  the  strain  energy  release  rates  Qq  and  corresponding  to  two  half 
plane  solutions  obtained  from  (54)  for  a  single  pressurized  crack  and  for  h  =  0  and 
h  =  X,.  respectively.  (Fig.  1). 

The  results  given  in  this  section  are  obtained  under  mode  I  loading  <7g  and  plane 
strain  conditions,  and  consist  of  largely  the  strain  energy  release  rate  §  which,  in  the 
case  of  the  interface  cracks,  is  accepted  to  be  the  main  crack  driving  force.  Figures  3 
and  4  show  some  sample  results  for  the  material  combinations  given  in  Table  2. 
Examples  F  through  H3  and  .1  correspond  to  structural  applications  involving  metals 


and  fiber  reinforced  composites  bonded  through  an  epoxy  layer.  In  I  the  adherents  are 
silicon  carbide  and  alumina  and  the  interfarial  region  is  amorphous  carbon.  In  examples 
Hj.  H2  and  H3  the  crack  is  between  aluminum  (med.  1)  and  epoxy  layer  (rued.  2)  and 
the  third  material  (med.  3)  is  a  fiber  reinforced  composite.  In  all  cases  the  variable  is 
h/2a.  Note  that  for  each  case  the  normalizing  stress  intensity  factor  Qy  is  given  in  Table 
2  and  is  calculated  from  (54)  (for  the  corresponding  material  pairs  1  and  3).  The  figures 
show  that  Q/Qq  approaches  1  for  (h/2a)— 0  and  for  (h/2a)— sc  (the  latter  limit  is 

also  given  in  Table  2).  Also,  since  the  stiffness  of  the  interfacial  layer  is  less  than  that  of 
the  medium  3.  Q  is  in  all  cases  a  monotonically  increasing  function  of  h/2a. 

Figures  5  and  6  show  the  strain  energy  release  rate  for  material  combinations 
K)....  K,  defined  in  Table  2.  .A.gain  Qq  corresponds  to  h  =  0.  Note  that  in  case  since 
the  stiffness  of  medium  2  is  greater  than  that  of  medium  3  we  have  >  0,^. 

The  results  for  the  second  group  of  examples  .A  through  E  described  in  Table  3 
are  shown  in  Figures  7-13.  The  parameter  e  shown  in  the  figures  corresponds  to  the 
ratios  £211 /E-j-  E222/E3  t)r  G212/G3  and  represents  the  structure  of  the  interfacial  zone. 
Here  again  Qg  is  the  strain  energy  release  rate  corresponding  to  h  =  0.  that  is.  to  the 
material  pair  1  and  3.  Note  that  the  crack  is  between  materials  1  and  2  and 
consequently,  as  seen  from  the  figures.  5  >  Qg  if  the  stiffness  of  the  interfacial  region  is 
less  than  that  of  medium  3.  i.e..  for  e  <  1  and  Q  <  Qg  for  e  >  1.  .A-lso.  in  all  cases  except  in 
case  B  g  is  bounded  by  Qg  and  Figure  13  shows  the  influence  of  the  modulus  ratio  e 
on  the  strain  energy  release  rate  for  h/2a  =  0.4.  For  =  1  materials  2  and  3  are  identical 
and  5  =  Qg.  The  figure  shows  that  as  expected  9/§g  is  a  monotonically  decreasing  function 
of  e.  The  exception  is  again  Case  B  in  which  the  modulus  E211  parallel  to  the  interface 
seems  to  have  very  little  influenc>='  on  Q.  This  result  may  also  be  observed  in  Figure  8. 


The  results  for  the  collinear  cracks  obtained  under  plane  strain  conditions  and 


uniform  tension  ctq  are  shown  in  Figures  14-20.  In  all  these  Figures  the  normalizing 
strain  energy  release  rate  Q"  is  obtained  for  a  single  crack  of  length  2a  and  for  the 
particular  material  or  material  combination  with  the  given  h/2a  ratios.  Figures  14  and 
15  show  the  results  for  an  isotropic  homogeneous  plane  (full  lines)  and  two  dissimilar 
isotropic  half  planes  (dashed  lines)  containing  two  interface  cracks  of  lengths  2a  and  2c. 
The  properties  of  the  dissimilar  half  planes  (E|  =  1  GPA.  E2  =  15  GPA.  i/,  =  =  0.1)  are 

selected  to  give  an  unusually  large  oscillation  parameter  (^- =  -0.1306).  In  both  cases 
the  limiting  values  of  the  normalized  strain  enei'gy  release  rates  are 

g,4/r  =  Qs/g’  =  1  .  =  Sd/Q'  =  .  for  d/a  . 

=  (c +  a)/a  .  for  d/a  =  0  .  (70) 

where  d  is  the  distance  between  the  inner  crack  tips  D  and  A  and  the  subscripts  C.  D. 
.A  and  B  refer  to  the  crack  tips  shown  in  the  figures,  c  =  a  in  the  example  shown  in  Fig. 
14  aiul  c  =  0.2a  in  Fig.  15.  For  the  dissimilar  materials  considered  we  have  g*  =  1.50363 
rr,j-a  I  GPA)''.  A  rather  interesting  and  very  useful  result  shown  by  Figures  14  and  15  is 
that  the  normalized  strain  energy  release  rates  for  the  collinear  cracks  in  bonded 
materials  may  be  approximated  by  that  obtained  from  a  homogeneous  isotropic  medium 
having  rhe  same  crack  geometry.  For  the  latter  very  often  it  is  possible  to  obtain  a 
closed  form  solution  (see.  for  example,  equations  (62)  and  (66)  which  were  used  to  plot 
the  solid  curves  in  Fig.  14). 

Figures  16-20  show  the  normalized  strain  energy  release  rates  in  bonded 
(Uthutropic  materials  (material  combination  G.  Table  2)  containing  two  collinear 
interface  cracks  with  c/a  =  1.  0.5.  0.2.  0.1  and  0.05.  respectively.  In  these  figures  solid 
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lines  again  correspond  to  the  homogeneous  isotropic  plane  with  the  same  crack 
geometry  as  the  bonded  planes.  The  round  dots  give  the  results  for  two  half  planes 
(h  =  30,  materials  a  and  d)  and  the  square  dots  for  two  half  planes  (materials  a  and  c) 
bonded  through  a  layer  (material  d)  with  h/2a  =  1.  It  may  be  seen  that  §/g‘  ratios  for 
the  two  dissimilar  half  planes  are  nearly  identical  to  the  isotropic  plane  values.  There  is. 
however,  some  discrepancy  between  the  isotropic  plane  and  the  interfacial  layer  r  nilts. 
Table  4  shows  the  limiting  values  of  Q/Q’  given  in  Figures  16-20  for  d/a  —0  and  d/a 
—  x:.  Note  that  for  h/2a=3c  the  limits  for  homogeneous  isotropic  and  bonded 
orthotropic  materials  are  identical,  whereas  for  h/2a  =  1  both  limit':  differ  from  the 
homogeneous  plane  results.  Despite  this  considering  the  addi'ionai  difficulties  involved 
in  solving  the  multiparameter  problem  of  collinear  cra.ks  in  bonded  materials,  the 
discrepancies  may  be  tolerable.  For  the  material  combination  G  and  for  various  h/2a 
ratios  the  value  of  Q’  are  given  in  T  'ble  5.  Some  additional  results  for  the  collinear 
cracks  for  a  fixed  d  (d/a  =  0.1)  and  v  rving  values  of  h/2a  .aay  be  found  in  Table  6. 

Some  sample  it  suits  fe:  the  stress  intensity  factors  along  with  the  corresponding 
strain  energy  rr  •  "  'ates  a  ^  t  in  Tables  7  and  8.  Equations  (48)  show  that  a.side 
from  the  material  constants  id  the  length  parameters,  the  stress  intensity  factors 
depend  on  the  im  ensionless  quantities  gj^.(  sf  1)  defined  in  (40).  In  the  examples 
considered  the  external  load  is  uniform  tension  at  infinity.  Thus,  the 

dimensionless  cprantities  F/cr^  given  in  Tables  7  and  8  to  represent  the  stress  intensity 
factors  at  the  crack  tips  aj,  bj,  a.^.  b^  (Fig.  1)  or  C.  D,  A.  B  (Fig.  17)  are  defined  by 

F(C)  =  3,,jr:7’ g.(-i).  F(D)=  g,(i). 


F(.4)  =  -,„jl-7^  g,(-l).  F(B)=  -7„>11-7^  g,(l). 

2.5 


(71) 


Note  that  7  is  dimensionless  and  721  has  the  dimension  of  stress.  The  material 
combination  G  is  again  used  in  the  examples  (see  Tables  1  and  2)  for  which  we  have 
721  =  1.8098  GPA.  7,2  =  1.7684  GPA,  a;  =  0. 071765.  From  (48)  and  (71)  one  may  als(; 
note  that  the  stress  intensity  factors  are  proportional  to  ^(b,^  -  a,j)/2  as  well  as  F.  Thus 
for  c/a  =  0.2  despite  the  seemingly  high  values  of  F(D)  gi\-en  in  Table  7.  from  a  fracrure 
propagation  view  point  by  far  the  worst  location  is  the  crack  tip  A.  This  may  l)e  seen 
from  the  values  of  the  strain  energy  release  rates  for  which  the  same  normalization 
constants,  o-^a  is  used  in  all  crack  tips.  In  the  case  of  collinear  cracks,  as  may  be  seen 
from  Figures  15  and  17-20,  this  conclusion  appears  to  be  quite  general  meaning  that,  for 
example,  under  cyclic  loading  fastest  subcriiical  crack  growth  would  take  place  at  the 
crack  tip  .A.  Also  it  appears  that  for  small  values  of  d/a  <  5b  <  §4-  '^’hereas  for 

greater  d/a  one  would  have  9c  <  <  §b  <  5.4-  This  may  have  bearing  on  the  crack 

uucleation,  growth  and  coalescence  ahead  of  a  main  crack  in  creep,  corrosion  as  well  as 
fatigue  crack  growth  studies. 
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Table  1.  The  elastic  constants  of 


GPA.) 


Mat.  a 

Mat.  b 

Ex 

153.09 

40.403 

E. 

40.405 

153.09 

E, 

22.754 

22.754 

Gxy 

29.304 

29.304 

1.551 

4.082 

Gx. 

4.082 

1.551 

Vxy 

1.834 

0.484 

t'xz 

0.195 

0.261 

^rjz 

0.261 

0.195 

Material  E  i/ 


d  (epoxy) 

3.1 

0.35 

e  (aluminum) 

69 

0.3 

f  (steel) 

200 

0.3 

g  (amorphous  carbon)  8.3 

0.32 

h  (silicon  carbide) 

207 

0.2 

i  (alumina) 

325 

0.3 

used  in  the  examples  (moduli  in  units  of 


Mat.  c 

Mat.  k 

Mat.  1 

30.6 

282.486 

225.734 

39.0 

225.734 

282.486 

6.4 

261.097 

261.097 

19.7 

131.06 

131.06 

4.5 

132.98 

108.11 

4.5 

108.11 

132.98 

0.351 

0.312 

0.312 

0.275 

0.244 

0.149 

0.275 

0.149 

0.244 

Material 

E 

u 

j  (high  modulus 

carbon)  380 

0.32 

m, 

2.5 

0.3 

m^ 

25 

0.3 

1113 

250 

0.3 

1114 

1000 

0.3 
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Table  2.  Material  Combinations  F-K4  used  in  the  examples,  o-  is  the  uniform  crark 
surface  pressure  and  2a  is  the  crack  length,  (see  Fig.  1  for  geometry  and  Table  1  for  the 
properties  of  materials  a-/) 


Mat. 

Med. 

Med. 

Med. 

g^/fT^a 

Fig. 

Comb. 

1 

2 

3 

(GPA)- 

'  (GP.A)- 

1 

F 

a 

d 

b 

0.0343 

0.4532 

13.2 

0.0718 

3 

G 

a 

d 

c 

0.0569 

0.4532 

7.966 

0.0718 

3 

H 

a 

d 

e 

0.0427 

0.4532 

10.62 

0.0718 

3 

Hi 

e 

d 

a 

0.0427 

0.4531 

10.62 

0.0671 

4 

H, 

e 

d 

b 

0.0320 

0.4531 

14.14 

0.0671 

4 

H3 

e 

d 

c 

0.0565 

0.4531 

8.027 

0.0671 

4 

I 

h 

g 

i 

0.0116 

0.1722 

14.88 

0.0843 

3 

.1 

e 

d 

f 

0.0273 

0.4531 

16.61 

0.0671 

3 

K, 

k 

m, 

1 

0.0109 

0.5501 

50.54 

0.0471 

5 

K, 

k 

1^2 

1 

0.0109 

0.0605 

5.60 

0.0741 

5 

K3 

k 

m3 

1 

0.0109 

0.0165 

1.06 

-0.0105 

6 

K, 

k 

1114 

1 

0.0109 

0.007 

0.645 

-0.0749 

6 
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Table  3.  Material  combinations  A-E  used  in  the  examples.  Referring  to  Fig.  1  and  Table 
1,  medium  3  is  epoxy  (d)  in  all  cases;  medium  1  is  amorphous  carl)on  (g)  in  A-E  and 
high  modulus  carbon  (j)  in  case  D,.  The  moduli  of  medium  2  are  selected  as  shown;  the 
Poissons  ratios  are  given  by  >^212^^^231  =  ^.23.2  =  0.35  =  The  external  load  is  unitorin 
crack  surface  pressure  <t  and  the  crack  length  is  2a. 


Case 

^211 

E., 

E222 

E.3 

G212 

9^/9o 

Fig. 

(GPA)-' 

(GPA)-’ 

A 

1 

1/4 

G3 

0.6139 

0.8622 

1.404 

0.0333 

7.  13 

1 

4 

G3 

0.6139 

0.3822 

0.623 

-0.0196 

B 

1/4 

1 

G3 

0.6139 

0.6056 

0.986 

0.0128 

8.  13 

4 

1 

G3 

0.6139 

0.5791 

0.943 

-0.0129 

c 

1/4 

1 

E211 

0.6139 

0.9745 

1.587 

0.0017 

9.  13 

2(1  +1/3) 

4 

1 

0.6139 

0.3995 

0.651 

-0.0115 

D 

1 

1/4 

E222 

0.6139 

1.4111 

2.299 

0.0208 

10.  13 

2(1+‘^3) 

1 

4 

0.6139 

0.2849 

0.464 

-0.0276 

E 

1 

1 

1/4  G3 

0.6139 

0.9792 

1.594 

0.0173 

12.  13 

1 

1 

4G3 

0.6139 

0.4576 

0.745 

0.0370 

D, 

1 

1/4 

E222 

0.4344 

1.2349 

2.842 

0.0396 

11.  13 

2(1 +^3) 

1 

4 

0.4344 

0.1167 

0.269 

0.0383 
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Table  4.  The  liiniriiia;  \'alu('s  uf  rlie  >riaiu  euei«y  n-lea.se  rates  slmwii  in  Fiiiures  lG-21). 

''‘)rrespouils  tu  material  (aniiljiiiatioii  G  «iveii  in  Table  2.  Luadimr  is  uiiiturm 
tension  rr|,.  ( Q*  ),„,=(  1 e  )iT,,‘a/S/i.  ( Q*  ().4o32()rr, pa  (GPA)''  anil 

( g* )  =  0.41  IbS  iT||‘a  (GPA)”'  an'  the  sin“;le  craek  vahu's.  tor  a  homoyen' ‘ohs 

meilium.  for  two  dissimilar  half  planes  (materials  a  and  dl.  and  tor  two  halt  plain"' 
bonded  throuyh  a  layer  (mat.  comb.  G.  h/2a  =  1).  respe<  ti\'ely. 

e/a 


d/a 

h/2a 

g/r 

1 

0.5 

0.2 

0.1 

0.05 

X 

- 

iQr/g*),,,, 

1 

0.5 

0.2 

0.1 

0.05 

X 

1 

1 

0.535 

0.219 

0.110 

0.055 

0 

- 

(58.1  vs*),.,,. 

o 

1.5 

1.2 

1.1 

1.05 

0 

1 

1.694 

1.372 

1.166 

1.0S3 

1.041 

X 

X 

(Qr/Q*).,-!/,. 

1 

0.5 

0.2 

0.1 

0.05 

0 

X 

l..r//i. 

o 

1.5 

1.2 

1.1 

1.05 

T.tbif  d.  \'alues  of  g*  for  the  material  eoinbination  G.  (Talde  3) 

h/2a  X  10  4  2  1  0.5  0 


0.453  0.4.53  0.450  0.440  0.411  0.34S  0.057 


;!i 


Q*^/a-||‘^a  (  GP.A. )  ' 


Table  6.  Strain  energy  release  rates  for  two  orthotropic  half  planes  bonded  through  a 
layer  of  thickness  h  (material  combination  G)  and  containing  two  collinear  cracks  of 
lengths  2a  and  2c.  d/a  =  0.1.  d  being  the  distance  between  the  inner  crack  tips  D  and  .A.. 
(Fig.  17) 

h/2a 


c/a 

OC 

10 

4 

2 

1 

0.5 

0 

1 

G4/cr“’a(GPA)-' 

1.489 

1.483 

1.457 

1.374 

1.142 

0.959 

0.187 

1 

G5/cT2a(GPA)-' 

0.605 

0.603 

0.592 

0.561 

0.484 

0.399 

0.076 

0.1 

G_j/<T'®a(GPA)'' 

0.533 

0.531 

0.528 

0.512 

0.456 

0.379 

0.067 

0.1 

Gg/<T^a(GPA)-' 

0.459 

0.458 

0.455 

0.444 

0.414 

0.349 

0.058 

0.1 

Go/cr^a(GPA)-’ 

0.209 

0.208 

0.206 

0.206 

0.191 

0.168 

0.026 

0.1 

G(;./(T‘*a(GPA)'’ 

0.135 

0.135 

0.134 

0.134 

0.126 

0.110 

0.017 

.T2 


Table  7.  Stress  intensity  factors  and  the  strain  energy  release  rates  in  ortlujtropic  half 
planes  bonded  through  an  isotropic  layer  and  containing  two  collinear  interface  cracks  of 
length  2a  and  2c.  c/a  =  0.2,  h/2a  =  l,  material  combination  G.  The  values  of  Q/c-^'a 
given  are  in  units  of  (GPA) 


d/a 

F(B)/<To 

F(A)/<T(j 

9,4/o-o‘a 

x; 

0.958 

-  0.172  i 

0.411 

0.958  -  0.172  1 

0.411 

4 

0.960 

-  0.172  i 

0.412 

0.961  -  0.172  , 

0.414 

2 

0.961 

-  0.172  ;■ 

0.413 

0.961  -  0.172  i 

0.414 

1 

0.961 

-  0.173  ; 

0.414 

0.962  -  0.173  / 

0.415 

0.5 

0.962 

-  0.175  i 

0.415 

0.985  -  0.176  i 

0.434 

0.1 

0.970 

-  0.178  / 

0.422 

1.105  -  0.203  i 

0.547 

d/a 

F(D)/(To 

9£)/'^o'^a 

F(C)/<ro 

cc 

0.958 

-  0.172  i 

0.082 

0.958  -  0.172  i 

0.082 

4 

1.035 

-  0.149  i 

0.095 

1.034  -  0.149  r 

0.095 

2 

1.072 

~  0.156  i 

0.102 

1.056  -  0.154  i 

0.100 

1 

1.179 

-  0.173  i 

0.123 

1.154  -  0.167  i 

0.118 

0.5 

1.332 

~  0.206  i 

0.158 

1.271  -  0.186  i 

0.143 

0.1 

2.053 

-  0.378  i 

0.379 

1.560  -  0.245  i 

0.216 

Table  8.  Stress  intensity  factors  and  strain  energy  release  rates  in  ortliotropic  half  planes 
bonded  through  an  isotropic  layer  and  containing  two  collinear  cracks  of  ecpial  length 
2a.  material  combination  G.  The  values  of  Q/o-^^a  given  are  in  units  of  (GPA)  ~  *. 


d/a 

=  0.1 

II 

h/2a 

F(B)/<To 

^a  F(A)/(ro 

F(B)/ffo^a 

1.170 -0.153i 

0.605 

1.834-0.252/ 

1.489 

1.012-0.145/ 

0.453 

10 

1.167-0.174/ 

0.603 

1.813-0.357/ 

1.483 

1.011-0.145/ 

0.453 

4 

1.154-0.178/ 

0.592 

1.796-0.363/ 

1.457 

1.008-0.147/ 

0.450 

2 

1.121-0.194/ 

0.561 

1.737-0.383/ 

1.372 

0.996-0.151/ 

0.440 

1 

1.036-0.203/ 

0.484 

1.566-0.422/ 

1.142 

0.958-0.170/ 

0.411 

0.5 

1.018-0.207/ 

0.399 

1.429-0.408/ 

0.959 

0.873-0.202/ 

0.348 

;54 


Av/o, 


F.g.  2  Normal  component  of  the  crack  opening  displacement  Av  =  u.,-  an 

.nterface  crack  in  two  bonded  ortho.ropic  half  spaces  for  various  values  of  v./r,  and 

are  respectively  the  normal  and  shear  components  of  crack  surface  tractions.  Medmm 
1  .s  materral  c,  medium  2  is  material  a  (Table  1),  c„  =  ar„/(v„v^.(l 

r\  n  — I- 


-  0.07505,  7?  =  =  1.057. 


0/0 


G/0 


E^=Z5 


Fig.  5  Strain  energy  release  rate  in  two  ortliotropic  half  planes  bonded  through  ai 
isotropic  layer  and  containing  a  pressurized  crack.  Material  combinations  Kj  and  K2 
(Table  2).  g/g^  =  1  for  h/2a  =  0. 


ic  layer  and  containing  a  pressurized  interface  crack.  Material  combination  A 
,  G212  =  G3,  e  =  E222/E'3- 


G/G 


G/0 


«  1/9 

Fig.  13  Dependence  of  Q  on  the  variable  e  for  a  single  interface  crack  in  bonded 
materials  described  in  Table  3.  h/2a  =  0.4;  e  =  i/E3  in  B  and  C,  e  =  E222/E3  D 


and  e  =  G212/G3  in  E. 


d/a 


Fig.  14  Strain  energy  release  rates  for  two  equal  collinear  uniformly  pressurized  cracks 
in  a  homogeneous  isotropic  medium  (full  lines)  and  in  two  bonded  dissimilar  isotropic 
half  spaces  (dashed  lines  E^/Ei  =  15,  j/j  =i/2  =  0.1.  w=  -0.1306).  In  each  case  Q'  is  the 
corresponding  single  crack  value  (obtained  from  d/a  =  oc).  g^/Q'-cc.  gg/g"-2  for  d/a-O 


Fig.  16  Strain  energy  relea.se  rates  for  two  collinear  cracks  of  equal  length  (c  =  a)  in  an 
isotropic  homogeneous  plane  (full  lines,  Q*  =  (1  +  Kjc-Qa/S/x).  in  two  bonded  half  planes 
(ineterials  a  and  d.  round  dots.  h/2a  =  oc,  Q' =  0.453(Ty^a  (GPA)-')  and  in  two  half 
planes  (materials  a  and  c)  Ijonded  through  a  hiyer  (material  d)(square  dots.  h/2a=l, 
Q'  =  0.411iTQ‘^a  (GPA)  ~  *).  e.xternal  load  is  uniform  tension  cr,2  =  o-q  at  infinity. 


d/a 


Fig.  17  Same  materials  and  geometry  a.s  in  Fig.  16,  c/a  =  0.5,  D  and  A  are  the  inner  and 
C  and  B  are  the  outer  crack  tips. 


Fig.  18  Same  as  Fig.  16,  c/a  =  0.2. 


Fig.  19  Same  as  Fig.  16,  c/a  =  0.1. 


Fig.  20  Same  as  Fig.  16,  c/  i  =  0.05. 


